Atomic calculations and search for variation of the fine structure constant in quasar 

absorption spectra 
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A brief review of the search for variation of the fine structure constant in quasar absorption 
spectra is presented. Special consideration is given to the role of atomic calculations in the analysis 
of the observed data. A range of methods which allow to perform calculations for atoms or ions with 
different electron structure and which cover practically all periodic table of elements is discussed. 
Critical compilation of the results of the calculations as well as a review of the most recent results 
of the analysis are presented. 

PACS numbers: 31.25.Eb, 31.25.Jf 



I. INTRODUCTION 

Theories unifying gravity with other interactions as 
well as many cosmological models allow for space-time 
variation of fundamental constants. Experimental search 
for the manifestation of this variation spans the whole 
lifetime of the Universe from Big Bang nucleosynthe- 
sis to the present-day very precise atomic clock exper- 
iments (see, e.g. reviews 0, An evidence that 
the fine-structure constant a (a — e 2 /hc) might be 
smaller in early universe has been found in the anal- 
ysis of quasar absorption spectra [3, 13, IE B S B @- 
The analyzed data which came from the Keck telescope 
in Hawaii included three independent samples contain- 
ing 143 absorption systems which spread over red shift 
range 0.2 < z < 4.2. The fit of the data gives Sa/a = 
(-0.543 ± 0.116) x 10~ 5 [13]. If one assumes the lin- 
ear dependence of a on time, the fit of the data gives 
d\na/dt = (6.40 ± 1.35) x 10 -16 per year (over time in- 
terval about 12 billion years). 

A very extensive search for possible systematic er- 
rors [ll[ has shown that known systematic effects can 
not explain the result. Although it is still not com- 
pletely excluded that the effect may be imitated by a 
large change of abundances of isotopes in last 10 billion 
years, it would need a very unlikely "conspiracy" between 
several elements. It had been checked that different iso- 
topic abundances for any single element can not imitate 
the observed effect. 

Recently two other groups of researchers fl2l [l3l . ll4T | 
applied the same method of the analysis to a different set 
of data obtained from the VLT telescope in Chile and re- 
ported no variation of a. There was an intensive debate 
in the literature about possible reasons for disagreement. 
It was argued in particular that at least part of the dis- 
agreement may be attributed to the spatial variation of 
a. This argument is based on the fact that Keck tele- 
scope is in the Northern hemisphere while VLT telescope 
is in the Southern hemisphere. 

The results of [13] were recently questioned in Refs. [H, 
[iH ] . Re-analysis of the data of Ref. [13] revealed flawed 
parameter estimation methods. The authors of [HI, [lfjj] 



claim that the same spectral data fitted more accurately 
give Sa/a = (— 0.64 ± 0.36) x 10~ 5 rather than Sa/a — 
(-0.06 ±0.06) x 10~ 5 as in Ref. 0. However, even this 
new result may require further revision. 

All the results discussed above were obtained with the 
use of the so called many-multiplet (MM) method which 
was proposed in Ref. 1171 . This method uses atomic cal- 
culations [H GJ, S3, ETUI [H, H, EH, SI, [23 to relate 
the change in the value of the fine structure constant to 
the change in the frequencies of atomic transitions. It is 
more than an order of magnitude more sensitive to the 
variation of the fine structure constant than the analy- 
sis of the fine structure intervals used for this purpose 
before. 

In present paper we review the methods of atomic cal- 
culations used in the search for the variation of the fine 
structure constant and present critical compilation of the 
most accurate results of the calculations and discuss some 
future directions. 



II. ATOMIC CALCULATIONS 

In atomic units a = 1/c, were c is speed of light, and 
a = corresponds to non-relativistic limit. Therefore, to 
reveal the dependence of atomic frequencies on a we need 
to perform relativistic calculations based on the Dirac 
equation. Doing this way we include leading relativis- 
tic corrections of the order (Za) 2 . The role of smaller 
corrections, such as Breit and quantum electrodynamic 
(QED) corrections will be discussed in section ITTEl 

It is convenient to present the dependence of atomic 
frequencies on the fine-structure constant a in the vicin- 
ity of its physical value ao in the form 



u>(x) = ujq + qx, 



(1) 



where ujq is the present laboratory value of the frequency 
and x = (a/ao) 2 — 1, q is the coefficient which is to be 
found from atomic calculations. Note that 



diu 
^ dx 



(2) 



x=0 
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TABLE I: Ab initio methods of atomic calculations depending 
on the number of valence electrons (N v ). 



Nv Method 



Accuracy 



1 MBPT + all-order sums 0.1 - 1% 

2-8 MBPT + Configuration interaction (CI+MBPT) 1-10% 
2-15 Configuration interaction (CI) 10-20% 



w 




To calculate this derivative numerically we use 
uj(+S) - u)(-6) 
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(3) 



and vary the value of a in the computer code. 

We use relativistic Hartree-Fock method as a starting 
point of all calculations. The radial equation for single- 
electron orbitals has the form (atomic units) 



dgv 
dr 



+ ^f v {r)- 2 + a 2 {e v -V HF ) g v (r)=0 



^Mr) + (e v - V HF )fv(r) = 0, 



(4) 



here k = (— 1) Z +J+ 1 / 2 (j + 1/2) , index v replaces the three- 
number set of the principal quantum number, and total 
and angular momentum: n,j, I. Vrf is the Hartree-Fock 
potential. 

Equations (@| with a = a VS + 1 are solved self- 
consistently for all core states to find Hartree-Fock po- 
tential of the atomic core. Then this potential is used 
to calculate a full set of single-electron orbitals for the 
states above the core. 

After single-electron states are calculated, the actual 
choice of the methods to calculate many-electron states 
of valence electrons depends on the number of valence 
electrons. Table Q] summarizes the methods used in the 
calculations. These methods will be discussed in more 
detail in following sections. 



A. Atoms with one external electron 

Atoms and ions of astrophysical interest which can be 
considered as systems with one external electron above 
closed shells in both ground and excited states include 
C IV, O VI, Na I, Mg II, Al III, Si IV, Ca II, Zn II and 
Ge II. All calculations for these atoms and ions were per- 
formed in the V N ~ l approximation in which Eq. ((4]) are 
solved for the closed-shell core with the external electron 
removed. The states of valence electron are calculated in 
the V N_1 potential of the frozen core. Correlations are 
also included using different numerical techniques. In 
most of the calculations (see, e.g. Ref. [3) the correla- 
tion potential method [281 ] was used. 

In this approach the correlation potential £, which de- 
scribes correlations between external and core electrons, 
is calculated in the lowest, second-order of the many- 
body perturbation theory (MBPT) using the B-spline 
basis set [29| . All four second-order diagrams for £ are 
presented on Fig. [TJ 





FIG. 1: Second-order diagrams for the matrix elements 
(ujsj 2 '' \w) of the single-electron correlation operator E^ 2 '. 



The £ operator is a non-local operator similar to the 
Hartree-Fock exchange potential. It is used in the equa- 
tions for the states of an external electron to calculate the 
so-called Brueckner orbitals (BO). The orbitals and cor- 
responding energies include all second-order correlations 
as well as higher-order correlations which correspond to 
(£) 2 , (£) 3 , etc. The accuracy for the energies in this 
approach is usually few per cent or better. 

Correlation potential method with the second-order £ 
does not include dominating higher-order correlations. 
The effect of the higher-order correlations on the q- 
factors for mono-valent atoms of astrophysical interest 
were studied in detail in Ref. [25| within the single-double 
coupled cluster approach and were found to be small. 



B. Atoms with open s and p shells 

Atoms or ions with two or more electrons on open 
s and/or p shells are next in the complexity of the 
calculations. Systems found in astrophysical observa- 
tions include C I, C II, C III, O I, O II, O III, 
O IV, Mg I and Al II. We perform calculations for 
such systems with the use of the configuration interac- 
tion method combined with the many-body perturba- 
tion theory (CI+MBPT) [30| . Configuration interaction 
(CI) technique is used for accurate treatment of the in- 
teraction between valence electrons while the MBPT is 
used to include correlations between core and valence 
electrons. Similar to the case of atoms with one ex- 
ternal electron, correlation operator £ is used to de- 
scribe the core- valence correlations. However, it now con- 
sists of at least two parts. Si describes correlations be- 
tween an external electron and the electrons in the core, 
while £2 describes correlation correction (screening) to 
the Coulomb interaction between two external electrons 
caused by the core electrons. The £1 operator is simi- 
lar to the £ operator used for atoms with one external 
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electron. The actual form for of the E operator depends 
on the choice of the potential in which core states are 
calculated. The simplest form is in the so-called V 
approximation [3lL |32| which is a generalization of the 
V N ~ X approximation used for atoms with one external 
electron. 

In this approximation the M external electrons are ex- 
cluded from the initial Hartree-Fock procedure for the 
core to make sure that the effective potential of the CI 
Hamiltonian for the valence states and the potential in 
which core electrons are calculated are the same. This 
is the key for the simplest form of the MBPT. This ap- 
proach gives good results for atoms with open s or p 
shells because electrons of these shells are localized on 
large distances and have very small effect on the wave 
functions of the core states (in spine of large effect on 
their energies, see Refs. [HI, for details). 

In the V N ~ M approximation the £1 operator is iden- 
tical to the £ operator used for atoms with one exter- 
nal electron. Corresponding diagrams are presented on 
Fig. Q] Diagrams for £2 are presented on Fig. [2] 

If the potential in which core states are calculated is 
different from the V N ~ M potential, then both Ei and E2 
include the so-called subtraction diagrams which account 
for this difference (see Ref. [30] for details). 

The effective CI Hamiltonian for the valence electrons 
has the form 



M 
i=l 



M 

£ 

i<3 



h E2ij 

Tii 



(5) 



where ho is a single-electron Hartree-Fock Hamiltonian 
which corresponds to Eq. ((4]). 
The CI equations 



H 01 ^(n,...,r M ) =E^{n,...,r M ) 



(6) 



are solved by matrix diagonalization using the the many- 
electron basis states constructed from the B-splines. 
Therefore, the B-spline basis set serves the dual purpose 
in the calculations: to calculate the correlation operator 
E and to solve the CI equations. 

The details of the calculations vary slightly from atom 
to atom. For example, E was not included for light atoms 
like carbon. Here relativistic corrections are small and 
high accuracy of the calculations is not needed. 



C. Atoms with open d or f shells 

Atoms with open d and / shells are most difficult for 
calculations. This is mostly due to large number of elec- 
trons which enter the CI calculations. This makes it prac- 
tically impossible to saturate the basis with an arbitrary 
choice of single-electron basis states. 

Also, accurate treatment of the core-valence correla- 
tions is much more difficult. First, the V approxi- 
mation is not valid because d and / electrons are inside 





FIG. 2: Second-order diagrams for double-electron correlation 

r (2) 

operator . 



the core and cannot be neglected in the initial HF cal- 
culations for the core. Therefore, subtraction diagrams 
must be included. Second, for atoms with more than 
two valence electrons, a new kind of E operator appear 
even in the lowest second-order of the MBPT, the three- 
electron operator £3 (see Ref. [30j for details). The most 
detailed study of an atom with open <i-shell, which used 
large basis, all three E-operators E1.E2 and E3 and Breit 
interaction were performed in Ref. [33j for the Fe II ion. 

We use a less sophisticated approach which still gives 
reasonably accurate results [H, [H, H3, HH, H3| • It is based 
on the CI technique with simplified treatment of the core- 
valence correlations (see below) and different choice of 
the basis in different calculations for the valent states. In 
our ear ly calculations [jj, M, HI H3 we used either a B- 
spline basis set or a set of orbitals which was constructed 
using a recurrent procedure suggested in Ref. [3~ij ] . In this 
procedure an virtual orbital is constructed from a lower 
one by multiplying it by distance r and orthogonalizing 
it to all lower orbitals (see also Ref. [30]|). 

In the most recent work [27j we used a more sophisti- 
cated approach in which a set of non-orthogonal Hartree- 
Fock valence states was constructed. The self-consistent 
HF procedure was done for each configuration of inter- 
est separately. Only one single-electron basis function in 
each partial wave is included in the CI for each config- 
uration. However, this functions represent good initial 
approximations since they come from the self-consistent 
HF calculations. 

To include core- valence correlations £1 is replaced by 
a correction to the local part of the HF potential of the 
atomic core 



SV = - 



2(r 4 + a 4 ) 



(7) 



Here a p is polarization of the core and a is a cut-off pa- 
rameter (we use a = as). The form of the SV is chosen 
to coincide with the standard polarization potential on 
large distances (— a p /2r 4 ). The a p is treated as a fitting 
parameter and values of a p for each configuration are cho- 
sen to reproduce their position in experimental spectrum. 
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The £2 and £3 operators are either neglected or simu- 
lated by introducing screening parameters for Coulomb 
integrals. It is usually assumed that the screening pa- 
rameters fk depend only on the multipolarity k of the 
Coulomb interaction and their values are chosen to have 
better agreement with the experimental energies. It is 
also assumed that all fitting parameters a p and fk do 
not depend on a. This is justified because the fitting 
represents only small correction to the energies. 
The effective CI Hamiltonian has the form 

M M 1 

i=l i<j y 

It can be obtained from ([5]) by dropping £2 and replacing 
£1 by SV [see Eq. ©]. 



where J, L, S are total and angular momentums and spin 
of the states. Eqs. (jlOllip can be used to find the mixing 
angle <j) which than is used to correct the calculated values 
of qi and qi. 

Note that it follows from ^ that the sum 51+52 — 
g( LS ) _\_g( LS ) anc j (j oes no -t; depend on <j). This can be used 
to check the accuracy of the two-level approximation. 

This scheme cannot be used if experimental values of 
g-factors are not known or if g\ + 52 7^ 9^^ + or if 
g( LS ) _ g^ LS \ If experimental g- factors are known but 

.9i + .92 7^ .9i LS ' + 9% then a multilevel consideration 
may help as explained in Ref. [2(J. Otherwise the only 
way to get reliable results is to fit the energies to experi- 
mental values to the accuracy sufficiently better than the 
small energy interval between strongly mixed states (see, 
e.g. [13). 



D. Level pseudo-crossing 



E. Breit and QED corrections 



Calculations may be complicated significantly by the 
phenomenon of level pseudo-crossing [2(| ■ In atoms 
with dense spectrum the states of the same total momen- 
tum J and parity and separated by small energy interval 
may be strongly mixed which in turn may lead to insta- 
bility of the calculations of the q coefficients. This may 
be considered as level pseudo-crossing in the vicinity of 
the physical value of a if energies of these two states are 
considered as functions of a 2 . In two- level approximation 
the q coefficients for each state are 



92 



Qi 



2 , (LS) . . 2 j (LS) 

qi = cos (j) q\ + sin (f> q 2 1 
{LS) +co S ^qi LS) 



(9) 



where <p is the mixing angle and q[ LS ^ and q^^ are the q 
coefficients for the states which correspond to the values 
of a far from the crossing. Very often atomic states in 
the absence of crossing are well described by the LS cou- 
pling scheme which is indicated by the (LS) superscripts 

in Eq. [5J If q[ LS ^ and q^^ are significantly different 
than the values of q± and q<i strongly depend on the mix- 
ing angle <j> which leads to instability of the results of 
calculations. 

The best way of dealing with this instability is to use 
experimental values of the Lande g-factors. In the two 
level approximation the (/-factors of two states are given 
by the formula very similar to Eq. @ : 



2 ± (LS) . ■ 2 j. (LS) 

gi = cos z (f) g\ +sm </> g 2 ; , 

• 2 , (LS) . 2 ± (LS) 

92 = sin 4> g\ +cos z </> g\ 



(10) 



with the same mixing angle cj> as in ([9]) and with g\ 
and g^^ given by 

(LS)_, , J(J + 1)-L(L + 1)-S(S + 1) 
9 + 2J(J + 1) ' 1 ' 



So far we have considered relativistic corrections to 
atomic energies which come from the Dirac equation. 
These corrections have an order of (Za) 2 . There are 
however other smaller corrections like Breit and quan- 
tum electrodynamic corrections. It is important to un- 
derstand the role these corrections may play in the search 
for variation of a. Breit corrections are proportional to 
a 2 but smaller than 2 power of nuclear charge Z. Its 
relative contribution is therefore larger for light atoms. 
In contrast, the QED corrections are proportional to a 3 
but stronger than Z 2 functions of Z. Their relative con- 
tribution is larger for heavy atoms. The role of the Breit 
interaction was studied in detail in Ref. [26] for mono 
and double-valent electron atoms and ions of astrophys- 
ical interest and in Ref. [HI for Fe II. The role of the 
QED correction has not been studied yet but expected 
to be small. 

The following form of the Breit operator is used in the 
relativistic calculations for many-electron atoms (atomic 
units) 



o.\ ■ a 2 + (&i ■ n)(&2 ■ n) 



(12) 



Here r = hr, r is the distance between electrons and 
a% is the a-matrix of the corresponding electrons. This 
is a low frequency limit of the relativistic correction to 
the Coulomb interaction between electrons. It contains 
magnetic interaction and retardation. 

It is important to include H B into Hartree-Fock Hamil- 
tonian to take into account the effect of Breit interaction 
on the self-consistent field (relaxation effect), for atoms 
with two or more external electrons Breit term is also in- 
cluded as a correction to the Coulomb interaction in the 
CI Hamiltonian. 

It turns out [26|, [33[ that Breit corrections to the q- 
values are relatively small and unlikely affect the analysis 
of quasar spectra in terms of variation of a. However, it 
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TABLE II: Energies of the transitions from the ground state 
of singie-valent atoms and ions of astrophysical interest and 
corresponding relativistic q coefficients (cm -1 ). 



Z Atom State E(expt) E(calc) 


q 


6 C IV 


2P1/2 


64484 


64504 


115(2 




2P3/2 


64591 


64636 


222(2 


11 Na I 


3P1/2 


16956 


16961 


45(0 




3P3/2 


16973 


16979 


62(0 




4P1/2 


30267 


30066 


57(1 




4p 3 /2 


30273 


30072 


51(1 


12 Mg II 


3^1/2 


35669 


35687 


121(1 




3P3/2 


35761 


35784 


212(1 




4pi/ 2 


80620 


80463 


161(1 




4p 3 /2 


80650 


80496 


192(1 


13 Al III 


3^1/2 


53684 


53723 


224(1 




3P3/2 


53917 


53970 


458(2 




4pi/ 2 


143633 


143538 


337(2 




4p 3 /2 


143714 


143623 


417(3 


14 Si IV 


3pi/ 2 


71290 


71352 


361(2 




3P3/2 


71750 


71836 


823(2 




4pi/ 2 


218267 


218226 


597(4 




4p 3 /2 


218429 


218397 


760(4 


20 Ca II 


4pi/ 2 


25192 


25086 


222(1 




4P3/2 


25414 


25315 


446(3 


30 Zn II 


4pi/ 2 


48481 


48721 


1541(7 




4f>3/2 


49355 


49606 2452(13 



is useful to include them for more accurate results. For 
example, as it was demonstrated in Ref. [26[ inclusion 
of Breit interaction significantly improves the agreement 
between theoretical and experimental fine structure mak- 
ing the result to be more accurate and therefore more 
reliable. 

Similar role is expected for the QED corrections. But 
it is the subject of future work. An adequate method to 
include leading QED corrections to the energies of many- 
electron atoms has been developed in Ref. [35| . 



III. RESULTS OF CALCULATIONS 

The results of the calculated g-values for lines which 
were observed in quasar absorption spectra are presented 
in Tables HH IIIII and IIVI Calculated and experimental 
energies are also presented to illustrate the accuracy of 
calculations. Experimental energies are taken from the 
NIST website |36j. Estimations of uncertainties for the 
values of q-factors are based on the sensitivity of the re- 
sults to variations of the calculation scheme and on the 
comparison of the experimental and theoretical energies 
and fine structure intervals. In cases when several cal- 
culations are available only latest most accurate results 
are presented. For the case of Fe II we also present the 
results of independent calculations by the St. Petersburg 
group [11]. 

The analysis of the values of the g-factors reveal an in- 
teresting picture. These values vary strongly from atom 
to atom and from one transition to another. The values 



can be very small or can be large positive or large nega- 
tive. This is due to the effect of several factors which can 
be qualitatively illustrated by the value of the relativistic 
energy shift of a single-electron state [l8| 



A„ = ^{Zaf 
v 



1 



i + i 



C(Z,j,l) 



(13) 



where E n is the energy of the single-electron state n, v is 
the effective principal quantum number (y = 1/ \J—2E n ), 
j is the total momentum of the state and constant C 
(C ~ 0.6) simulates the effect of exchange interaction 
and other many-body effects. The value of the g-f actor 
for a particular atomic transition can be roughly esti- 
mated assuming that the transition involves the change 
of a state of just one electron. Then 



(14) 



where state n' is above state n. 

The first and most obvious feature of the g-factors is its 
Z 2 dependence on the nuclear charge Z. However, it also 
depends, e. g. on the ionization degree. This is because 
the absolute value of the single electron energy E n (fT3"|) 
is larger for ions. Also, dependence of the energy shift 
(fT"3")) on the total electron momentum j suggests that the 
value of q is likely to by positive for the s — p transitions 
and negative for the p — d transitions. 

All these features of the g-factor values are reproduced 
in the accurate numerical calculations (see Tables IIIIIIII 
and llVp . This complex dependence of atomic frequencies 
on a is very important for study of possible systematics. 
It is very unlikely that any other unknown effect exhibits 
exactly the same features. Therefore, if the results of 
the analysis based on different transitions give the same 
result for variation of a then the mimic of the effect by 
any systematic is very improbable. 

Theoretic uncertainty for the values of q varies sig- 
nificantly from atom to atom. The most accurate re- 
sults are for atoms which can be considered as atoms 
with single valence electron above closed shells in both 
ground and excited states. The results for such atoms 
are presented in Table [TXJ They were obtained using 
second-order MBPT [26[ as well as the single-double 
coupled cluster approximation combined with the third- 
order MBPT [25|. Small uncertainty is due to excellent 
agreement between different approaches and between ex- 
perimental and theoretical data for the energies and fine 
structure intervals. Breit interaction was also included 
in Ref. [26| , Although Breit contributions are very small 
for the purposes of the analysis they are larger than the 
uncertainty of the calculation of the correlations. It was 
demonstrated in Ref. (2|| that inclusion of Breit con- 
tributions improve significantly theoretical fine structure 
brining it to almost perfect agreement with the experi- 
ment. 

Table [TTT1 presents results for Fe I and Fe II. Iron is one 
of the most important elements used in the analysis. It 
is used more often than most of other elements and it 
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TABLE III: Energies of the transitions from the ground state of Fe I and Fe II and corresponding relativistic q coefficients 
(cm- 1 ). 

Atom Ground state Upper state Energy q 



or ion Coring. Term 


J Config. 


Term 


J 


Expt.f36] 


Calc. 


[27] 


Tool 

[33J 


T T T O it) A A b 

Fe I 3d 4s D 


A O Tt> A A 

4 3d 4s4p 


b D° 


4 


25899 


26428 


999(300) 






3d 6 4s4p 


5po 


5 


26874 


27432 


880(260) 






3d 6 4s4p 


P 


3 


29056 


29340 


859(260) 






o to a 2, a 

3d 4s 4p 


5 D° 


4 


33095 


32680 


2494(750) 






O TO A 2 A 

3d 4s 4p 


5 D° 


3 


33507 


33134 


3019(900) 






O TO A 2 a 

3d As 4p 


5po 


5 


33695 


32522 


2672(800) 






3d 4s 4p 


JJ 


4 


39625 


OAT A A 

39544 


1680(500) 






3d 5 4s 2 4p 


5po 


5 


40257 


40194 


1042(300) 




Fe II 3d 6 4s 6 D 


9/2 3d 6 4p 


6 D° 


9/2 


38458 


38352 


1330(150) 


1410(60) 




3d 6 4p 


6 D° 


7/2 


38660 


38554 


1490(150) 


1540(40) 




3d 6 4p 


6po 


If/2 


41968 


41864 


1460(150) 


1550(60) 




3d 6 4p 


6po 


9/2 


42II4 


42012 


1590(150) 


1660(60) 




3d 6 4p 


6po 


7/2 


42658 


42715 


1210(150) 


1540(400) 




3d 6 4p 


4po 


7/2 


62065 


65528 


1100(300) 


1560(500) 




3d 5 4s4p 


6po 


7/2 


62171 


65750 


-1300(300) 


-1030(300) 



is also most studied theoretically. First analysis of the 
quasar absorption spectra compared shift of frequencies 
of Fe II to those of Mg I and Mg II. The values of the q 
coefficients for magnesium are small compared to those 
for Fe II. Therefore, one can say they were used as "an- 
chors" against which the shift of Fe II lines was measured. 
When it was realized that Fe II has also large negative 
shifters (see last line of Table llllj) it was suggested that 
Fe II alone can be used in the analysis by comparing lines 
with positive and negative values of q. This may help to 
eliminate certain types of systematic errors [33[ . 

Calculations for Fe II were carried out by means of the 
CI method in Refs. [H, H(|. The more accurate results 
from the later work are included in the Table. The Table 
also presents the results of independent calculations by 
the St. Petersburg group This is the most detailed 
and accurate calculations which include core- valence cor- 
relations and Breit corrections. The results of both cal- 
culations agree within the declared accuracy. Accuracy is 
high for lower states but deteriorates significantly higher 
in the spectrum due to increasing configuration mixing 
and difficulty in the achieving of basis saturation for as 
many as seven external electrons. 

Lines of neutral iron has not been used in the analysis 
yet. However, the lines presented in Table IIIII are ob- 
served in the quasar absorption spectra and are planed 
to be included into analysis. Therefore we have recently 
calculated the q coefficients for them 27]. It turned out 
that the values of q are large and for some lines they are 
even larger than for Fe II. This is due to strong configura- 
tion mixing which makes the transitions to be effectively 
two-electron transitions (see Ref. [27| for details). This 
makes Fe I a good candidate for the inclusion into the 
analysis. 

The results for other atoms are presented in Table IIVI 
The calculations were carried out in works of Refs. [la. 
0, H SI [H Hi. See these works for more detailed 



information. Here we present only the most commonly 
used lines and the latest more accurate calculations. 



IV. FUTURE DIRECTIONS 

The analysis of quasar absorption spectra has already 
produced very important results putting strong con- 
strains on possible space-time variation of the fine struc- 
ture constants a and hinting that a might be smaller 
in early epoch. However the results are controversial and 
more work is needed to prove or dismiss the contradicting 
claims. Resolving the disagreement between the analysis 
of the Keck and VLT data seems to be the most impor- 
tant thing at the moment. This should involve careful 
cross re-analysis of the data by independent groups of 
experts. 

It is also important to include more data into the anal- 
ysis to improve statistics significantly. For example, hav- 
ing sufficient statistics for different lines with different 
shifts would allow to compare the results which come 
from the analysis of each line separately and therefore ex- 
clude many possible systematics. Also, improving statis- 
tics may allow to study variation of the fine structure 
constant as a function of the red-shift parameter z or as 
a function of the position in the sky (space variation) . 

There are probably many more lines observed in quasar 
absorption spectra but still not used in the analysis. In- 
clusion of this data may also add important information 
about a variation. For example, calculations for lines 
of Fe I which are available for the analysis reveal that 
the frequencies shifts for these lines due to change of 
a are large and significantly different for different lines. 
Therefore, if corresponding frequency shifts are observed 
it would be hard to attribute them to anything but vari- 
ation of a. 

Another example is the inclusion of the weak Ml or E2 
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TABLE IV: Energies of the transitions from the ground state of atoms and ions of astrophysical interest and corresponding 
relativistic q coefficients (cm" 1 ). 



Z Atom 


Ground 


Upper 


Ener 


gy 


q 


or ion 


state 


state 


Expt.[36] 
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1 nni 7n 
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loU^DU ) 


p. r 1 tt 

O XX 


ZS zp 


2po 
r l/2 


zs zp 


2po 
r 3/2 


00 


7 1 


Oo(,lJ 






zszp 


D 5/2 


*7A AQA 


I doUd 


1 '70/ r >A\ 

1 <y(zu ) 








zszp 


U 3/2 


7 jnQQ 

1 4yoo 


( DOUo 


1 /0(zU j 








9q9<T1 2 

z^zp 


°l/2 


Afi /! A/! 


A7AAQ 

y ( yyo 


1/?1 / QA\ 
lDl(jU j 


a ( ' TTT 
U 111 


2s 2 


^0 


2s2p 


lpo 
r 1 


lUzOOZ 


1 aqoe;p; 


10o(^l j 


7 N V 

1 IN V 


2 s 




zp 


2po 
*3/2 


80799 
ou f 


SI RD7 


409^0^ 


S D TT 

own 


zs zp 


°3/2 


zszp 


r 5/2 


1 1 Q87^ 
1 lyo 1 


1 99R90 
rzzuzu 








ZsZp 


4p 

T3/2 


IzUUUU 


Izz / Do 


4oy(t>u j 








ZSZp 


" 1/2 


1 OflAC Q 

IzUUoa 


lzzo4o 


( 4(oU J 


Q ( \ TTT 

o U 111 


Zi zp 




9q9« 3 
z^zp 




IzUUoo 


izizyy 


( zo(t>U ) 








2s2p 3 


3po 
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/ ZO^OU ) 


s n tv 
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8 D VT 


2s 


2q 
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yuo 1 
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2d 


2po 
r 3/2 


yoyuo 


y < uy 1 


y lo^ou ) 


12 Mg I 


3s 2 




3s3p 


lpo 
A 1 


OULM 1 


oouou 










3s4p 


lpo 


49347 


49277 


soni 

olhi; 


13 Al II 


3s 2 




3s3p 
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59852 


59800 


270(1) 


14 Si II 


3s 2 3p 


2po 

1/ J 


3s3p 2 




55309 


54655 


520(30) 






3s 2 4s 


- ' 


65500 


54675 


50(30) 


22 Ti II 


3d 2 4s 


4 F3/2 
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396(50) 
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29401 
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4 T?° 
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0/2, 
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31143 
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25 Mn II 3d 5 4s 


7 s 3 
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*3 
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M 
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1276(150) 








3d 4 4s4p 


7 P§ 
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83363 
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3d 4 4s4p 


7po 
r 3 
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3d 4 4s4p 


7po 
M 
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28 Ni II 


3d 9 


2 D 5/2 


3d 8 4p 


2™ 
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57080 
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32 Ge II 


4s 2 4p 


3po 
r l/2 


4s 2 5s 


Sl/2 


62402 


62870 


-630(40) 



transitions between the states of the same fine structure 
multiplet (37I ]. It was argued that the anomalies in the 
fine structure intervals lead to enhancement of the sensi- 
tivity of the transition frequencies to variation of the fine 
structure constant. This idea is similar to what what 
suggested before for laboratory experiments [38[. 

It is also important to include molecular lines into anal- 
ysis. This would allow to study variation of the ratio of 



the electron to proton mass Most of grand unifica- 
tion models suggest that the variation of mass ratio and 
a are related and mass ratio is changing faster than a. 
Therefore it might be easier to find manifestations of vari- 
ation of mass ratio. There are already published results 
of such analysis performed by different groups [1^, H3] 
but as for the case of a the results are controversial but 
the other way around. Those who claim variation of a see 
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no variation of /x and vise versa. Here again more study 
is needed and there are some new interesting suggestions 
(see, e.g. [Hj]). 

V. CONCLUSION 

The many-multiplet method which is based on the 
analysis of the frequencies of strong electric dipole tran- 
sitions in atoms and ions found in gas clouds intersecting 
the sight line from Earth to distant quasars prove to be a 
useful tool for the search of the variation of the fine struc- 
ture constant in quasar absorption spectra. The method 
relies on atomic calculations to reveal the dependence of 
atomic frequencies on the fine structure constant. Criti- 
cal compilation of the all relevant calculations performed 
by our and some other groups and presented in this paper 
can serve as a reference point for future analysis. The re- 
sults of the analysis so far is controversial. The analysis 



of the data from Keck telescope in Hawaii indicate that a 
might be smaller in early Universe while similar analysis 
of the data from the VLT telescope in Chile performed 
by different groups of researchers gives null result. Both 
analysis use the same MM method and the same atomic 
calculations. This means the the reason for disagree- 
ment is probably not relevant to atomic calculations and 
is rather in the data or the analysis. Recent re-analysis 
of the VLT data performed by the authors of the anal- 
ysis of the Keck data pointed to some serious problems 
in the method used by the other group. Further revision 
of the VLT data is needed to resolve all problems in the 
analysis and disagreement in the results. 
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